MULTIDIMENSIONAL YAMADA-WATANABE THEOREM 
AND ITS APPLICATIONS 
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Abstract. Multidimensional and matrix versions of the Yamada-Watanabe theorem 
are proved. They are applied to particle systems of squared Bessel processes and to 
matrix analogues of squared Bessel processes: Wishart and Jacobi matrix processes. 



1. Introduction 

In this note we prove a multidimensional and matrix analogues of the celebrated 
Yamada-Watanabe theorem, ensuring the existence and uniqueness of strong solutions of 
one-dimensional stochastic differential equations (SDKs) with a Holder coefficient in the 
ltd integral part. In Section [2] we prove a multidimensional Yamada-Watanabe theorem 
(Theorem [2]). In Section [3] we show the existence and uniqueness of a strong solution of 
the matrix SDKs 

dXt = g{Xt)dBMXt) + h{X,)dB^g{Xt) + b{Xt)dt 

where g, h,b : H ^ H are such that g ® h is 1/2-Holder continuous symmetrized g"^ ® h'^ 
and b are Lipschitz continuous and Bt is a Brownian p x p matrix, and establish in this 
way a matrix Yamada-Watanabe theorem, see Theorem [5l Section |4] contains interesting 
applications. We apply Theorems [2] and |5] to 

(i) noncoUiding particle systems of squared Bessel processes which are intensely studied 
in recent years in statistical and mathematical physics (Katori, Tanemura [131 113] )• 

(ii) the systems of SDKs for the eigenvalues of Wishart and Jacobi matrix processes, 
as well as to the /3- Wishart and /3- Jacobi processes. We note the importance of the 
/9- Wishart systems in statistical physics: they are statistical mechanics models of "log- 
gases", see the recent book of Forrester [TT] . 

Surprisingly, the existence of strong solutions of SDKs for such "Holder" non-colliding 
particle systems was not established in general; only some cases of (ii) were treated by 
Demni [3, [8] . In Sections 14.11 and 14.41 we prove the existence and uniqueness of a strong 
solution to all these systems of SDKs. 

The matrix Yamada-Watanabe theorem is applied in Sections 14.21 and 14.31 to some ma- 
trix valued squared Bessel type processes. We improve the known results of Bru [21 13], 
Mayerhofer et al. [T7] and Doumerc |1U] on the existence and uniqueness of strong solu- 
tions to Wishart and Jacobi matrix SDKs. We extend them to the whole range of the drift 
parameter b{Xt). In the Wishart case we contribute in this way to realization of a pro- 
gramme started by Donati-Martin, Doumerc, Matsumoto, Yor [9] , claiming that Wishart 
processes have similar properties as classical 1-dimensional squared Bessel processes. 
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2. A MULTIDIMENSIONAL YAMADA-WATANABE THEOREM 

Let us recall the classical Yamada-Watanabe theorem, see e.g. [12], p. 168 and 

Theorem 1. Let B(t) be a Brownian motion on R. Consider the SDE 

dX{t) = a{X{t))dB{t) + b{X{t))dt. 

If |cr(x) — cr(y)p < p{\x — y\) for a strictly increasing function p on with p(0) = 
and J^^ p~^{x)dx = oo, and b is Lipschitz continuous, then the pathwise uniqueness of 
solutions holds; consequently the equation has a unique strong solution. 

No multidimensional versions of the Yamada-Watanabe theorem seem to be known, 
even if their need is great (cf. P- 738). We propose a useful generalization, however we 
stress the fact that the Holder continuous functions (jj appearing in the following system 
of SDEs are one- dimensional. The proof is based on the approach presented in Revuz,Yor 
|20) . in particular on the results of Le Gall [16]. By || ■ || we mean the Euclidean norm 
II ■ II 2 on R"^. 

Theorem 2. Let p,q,r G N and the functions bi : R^ — )■ R, i = 1, . . . ,p and Ck,dj : 
^p+r Yi, k = p + 1, . . . ,p + q, j = p + 1, . . . ,p + r, be bounded real-valued and 
continuous, satisfying the following Lipschitz conditions 

\hiyi) - bi{y2)\ < A\\yi - yall, i = l,...,p, 

\ckiyi,zi) - Cfc(?/2, 2^2)1 < A\\{yi,zi) - (1/2,22)11, k = p+l,...,p + q, 

\dj{yi,zi) - dj{y2,Z2)\ < A\\{yi, Zi) - (^2,22)11, j = p + 1, . . . ,p + r, 

for every yi, y2 G R^ and Zi, Z2 G R*". Moreover, let : H. H, i = 1, . . . ,p, be a set of 
bounded Borel functions such that 

\ai{x) - ai{y)\^ < p^{\x - y\), x,yGR, 

where pi : (0, 00) — )■ (0, cxd) are Borel functions such that Jq+ Pi^^{x)dx = 00. Then the 
pathwise uniqueness holds for the following system of stochastic differential equations 

dYi = (Ti{Yi)dB, + bi{Y)dt, z = l,...,p, (2.1) 

p+g 

dZj= Ck{Y,Z)dBk + dj{Y,Z)dt, J =p+l,...,p + r. (2.2) 

k=p+l 

Proof. Let (Y, Z) and {Y, Z) be two solutions with respect to the same Brownian motion 
B = {Bi)i<p+q such that Y{0) = Y{0) and Z{0) = Z{0) a.s. For every i = 1, ... ,p we 
have 

Yi{t) - Y,{t) = [\a,{Y,{s)) - a,{Y,{s)))dBi{s) + f\b,{Y{s)) - b,{Y{s)))ds. (2.3) 
Jo Jo 

Then we get 

„ p,(Y,{s) - y.M) ^ ' h Pi(Y,(s) - Y,(s)) W,)>yM) - 

Thus, applying Lemma 3.3 from [20] p. 389 , we get that the local time of Yi — at 
vanishes identically. Consequently, by Tanaka's formula we get 

\Y,{t)-Y,{t)\= I sgn(F,(s)-F,(s))rf(r,(s)-F,(s)) + L°(F,-r,) 







sgn{Y,{s) -Y,{s))d{Y,{s) -Y,{s)). 
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Since ai is bounded, using (12. 3p . we state that 

\Y,{t)-Y,{t)\ - [\gn{Y,{s) -Y,{s)){b,{Y{s)) -b,{Y{s)))ds 
Jo 

is a martingale vanishing at 0. This together with the Lipschitz conditions satisfied by 6, 
give 

E\Y,{t)-Y,{t)\<A [ E\\Y{s)-Y{s)\\ds. 
Jo 

Summing up the above-given inequahties we arrive at 

E\\Y{t)-Y{t)\\<C fE\\Y{s)-Y{s)\\ds 

Jo 

and Gronwall's lemma shows that Y{t) = Y{t) for every t > a.s. 

Using in a standard way the properties of the Ito integral and the Schwarz inequality, 
similarly as in [12], p. 165, we get that for every t G [0,T] 

E\Z,{t)-Z,{t)\^ <C 5^ E(/ {ck{Y{s),Z{s))-ci,{Y{s),Z{s))dBk{s)f 

k=p+i 

+ CE{ [\d,{Y{s), Z{s)) - d,{Y{s), ~Z{s))dsf 
Jo 

p+1 pt 

<C Y.E {ck{Y{s),Z{s))-Ck{Y{s),Z{s))fds 
k=p+i 

+ CTE [ {d,{Y{s),Z{s))-dj{Y{s),Z{s))yds 
Jo 

< CA\q + T)E !\\\Y{s) - Y{sW + \\Z{s) - Zis^ds 
Jo 

Thus, using the previously proved fact that Y = Y a,.&. we get that 

E\\Z{t) - Z{t)\\' < CA^{q + T)r f E\\Z{s) - Z{s)\\^ds. 

Jo 

One more application of the Gronwall's lemma ends the proof. □ 

3. Matrix stochastic differential equations 

Consider the space Sp of symmetric p x p real matrices. Recall that if : R — )■ R and 
X ^ Sp then g{X) = Hg{A)H^ , where X = HAH^ is a diagonalization of X, with H an 
orthonormal matrix and A a diagonal one. Denote by Bt a Brownian pxp matrix. Let Xt 
be a stochastic process with values in Sp such that Xq G Sp, the set of symmetric matrices 
with p different eigenvalues. Let = diag{Xi{t)) be the diagonal matrix of eigenvalues 
of Xt ordered increasingly: Xi{t) < A2(t) < . . . < Ap(t) and Ht an orthonormal matrix of 
eigenvactors of Xt. Matrices A and H may be chosen as smooth functions of X until the 
first collision time r = inf{t : Xi{t) = Xj{t) for some i ^ j}, ci. |19) . 

3.1. Eigenvalues and eigenvectors of Xf In order to prove a matrix Yamada- 
Watanabe theorem we need to consider the SDKs satisfied by the processes of eigenvalues 
and eigenvectors of Xt. 
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Theorem 3. Suppose that an Sp-valued stochastic process Xt satisfies the following ma- 
trix stochastic differential equation 

dXt = g{Xt)dBth{Xt) + h{Xt)dBjg{Xt) + h{Xt)dt (3.1) 

where gf, /i, 6 : R — )■ R, and Xq G Sp. 

Let G{x,y) = g'^{x)h'^{y) + g'^{y)h'^{x) . Then, for t < t the eigenvalues process At and 
the eigenvectors process Ht verify the following stochastic differential equations 

d\ = 2g{\)h{Xi)du, + (h{\) + T^T^ ) (3.2) 

dh, = j: ft, v|SlMi^,^ _ J2 §^,dt (3.3) 

where {ui)i and {Pkj)k<j o-f^ independent Brownian motions and = Pkj- 

Proof. The proof generahzes ideas of Bru [I] in the case of Wishart processes. Following 
|10| in the case of matrix Jacobi processes, it is handy to use the Stratonovich differential 
notation X o dY = XdY + ^dXdV. We then write the Ito product formula 

d{XY) = dX oY + X odY. 

We also have dX o (YZ) = {dXoY)oZ and (X o dY f = dY^ o X^. 
Define A, a stochastic logarithm of H, by 

dA = R-^ odH = H^ o dH. 

Observe that by Ito formula applied to H^H = I, the matrix A is skew-symmetric. By 
Ito formula applied to A = H^XH, setting dN = o dX o iJ, we get 

dk = dN + ko dA - dAo K. 

The process A o dA — dA o A is zero on the diagonal. Consequently dXi = dNa and 
= dNij + {Xi — Xj) o dAij, when i ^ j. Thus 

dA,j = ^ o dNij, i 7^ J- (3.4) 

Aj — Ai 

For further computations we need the quadratic variation {Xst^Xs/f/) which is easily 
computed from (13. ip : 

dX^tdX^^t, = g\X)s,,h\X)u' + g\X)s,h\X),,t + g\X),th\X),e + g\X)u'h\X),,.. 
The martingale part of dN equals the martingale part of H^dX H and by the last formula 

dNijdNkm = g\A)ikh\A)jm + g\A)imh^{A)jk + g^{A)Jkh\A)^n, + g\A)jmh\A)ik. (3.5) 
From (13. sp it follows that 

dN.dNjj = A5ijg\X,)h\X,)dt. (3.6) 
Now we compute the finite variation part dF of dN 

dF = H^b{X)Hdt + ^{dH^dX H + H^dXdH) 

= b{A)dt + ^ {{dH^ H) {H^dX H) + {H^dX H) {H^dH)) 

= h{A)dt + -{dNdA + {dNdAf). 
2 
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Using ( I3.4p and f l3.5p we find, writing G{x,y) = g'^{x)h^{y) + g'^{y)h'^{x) 



{dNdA)ij = ^ dNikdAkj = Sij ^ ^ _^ . 

k^j kj^i * 

It follows that the matrix dNdA is diagonal, so also dF is diagonal, 

dFu = h{K)dt + V y-^ dt 

and by (13 ■4p . A is a martingale. Finally, using (13. 6 p and the last formula, there exist 
independent Brownian motions z/j, z = 1, . . . , m, such that (13.21) holds. 
In order to find SDEs for iff, we deduce from the definition of dA that 

dH = HodA = HdA + -dHdA = HdA + -HdAdA. 

2 2 

By ([S3]) we find dN^dNij = g'^{Xi)h'^{Xj) + g^{Xj)h'^{Xi) and dN^dNkm = when the 
ordered pairs « < j and k < m are different. We infer from (13.40 that 

c^A.- = (3.7) 

where the Brownian motions (Aj)i<j are independent and = Moreover, when 
k < m, we have dXidakm = dNadNkm/ {Xm — Xk) = by (13. 5p . so the Brownian motions 
Wij)i<j ctnd (i/i)j are independent. From (13.70 we deduce that the matrix dAdA is diagonal 
and 

{dAdAh = -J2 dAkdAk = -J2 ^T%t4- 

k^t k^i 

Now we can compute dH = HdA + ^HdAdA and prove (13. 3p . □ 

3.2. Complex case. In this subsection we study the eigenvalues process for a process 
Xt with values in the space Tip of Hermitian p x p matrices. 

Theorem 4. Let Wt be a complex p x p Brownian matrix (i.e. Wt = B] + iBf where Bj 
and Bf are two independent real Brownian squared matrices). 

Suppose that an Hp-valued stochastic process Xt satisfies the following matrix stochastic 
differential equation 

dXt = giXt)dWth{Xt) + hiXt)dW:g{Xt) + 6(Xf)rft, (3.8) 

where g, h,b : H. ^ H, and Xq G Hp. 

Let G{x,y) = g^{x)h'^{y) + g'^{y)h'^{x) . Then, for t < t the eigenvalues process At 
verifies the following system of stochastic differential equations 

d\ = 2g{X,)h{X;)du;+ U(A,) + 2^ ^^^^'^'M rft, (3.9) 
where are independent Brownian motions. 

Proof. We will need the following formula for the quadratic variation [Xst^iXs't') which 
is computed from (13. 8p . using the fact that for a complex Brownian motion Wt, the 
quadratic variation processes satisfy d(W., W) = and d{W, W) = 2dt. 

dXstdX,.t' = 2 {g\X)st'h\X)s^t + g\X),,th\X)st') . (3.10) 
Define A, a stochastic logarithm of H , by 

dA = H~^ odH = H* o dH. 
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By Ito formula applied to H*H = I, the matrix A is skew-Hermitian. In particular, the 
terms of diag{A) are purely imaginary (recall that in the real case they were 0). By 
Ito formula applied to A = H*XH, we get, setting dN = H* o dX o H 

dA = dN + AodA-dAo A. 

We have 

dN = H*dXH + ^{dH*dX H + H*dXdH) 

so the process takes values in Tip. In particular its diagonal entries are real. The 
process A o dA — dA o A is zero on the diagonal, so dXi = dNa. Moreover, when i ^ j, we 
have = dNij + (Aj — Xj) o dAij and 

dAi, = ^ o dNij 1 i^]- (3.11) 

Aj — Ai 

The martingale part of dN equals the martingale part of H*dX H and by formula fl3.10p 
we obtain 

dN.^dNkm = 2 [g\A),rah\A)jk + g\A)^kh\A),„) • (3.12) 
From it follows that 

dN.dNjj = 4:5ijg^{Xi)h^{Xi)dt. (3.13) 
Now we compute the finite variation part dF of dN 

dF = H*b{X)Hdt+^{dH*dX H + H*dXdH) 

= b{A)dt + ^ {{dH* H){H*dX H) + {H*dX H){H*dH)) 

= h{A)dt + ]^{dNdA + {dNdA)*). 
Recall that G{x,y) = g^{x)h^{y) + g^{y)h'^{x). We get 

(dNdA),, = ^ dN^kdAk, = 25,, ^^^Z^^ + dN^jdA,,. 

k k^i ' ^ 

When i = j, the term dNa is real and dAa G zR. It follows that 

dF,, = b{X.)dt + 2Y2^^^^dt. 

k^i ^'-^'^ 

Finally, using fl3.13p and the last formula, there exist independent Brownian motions i/j, 
i = 1, . . . , m, such that (13. 9p holds. □ 

The Theorem H] may be applied in a special case g{x) = y^, h{x) = 1 and b{x) = 
26 > 0, when the equation (13.80 is the SDK for the complex Wishart process, called also 
a Laguerre process. This process and its eigenvalues were studied by Konig-O'Connell 
|15| and Demni [3]. 

3.3. Collision time. In this subsection we show that under some mild conditions on the 
functions g and h in the SDK (13. ip . the eigenvalues of the process Xf never collide. 

Proposition 1. Let A = (Aj)j=i,,,p be a process starting from Ai(0) < . . . < Ap(0) and 
satisfying 113. 2\} with functions b, g, h : H H. such that b, g"^, h"^ are Lipschitz continuous 
and g'^h'^ is convex or in class C^'^ . Then the first collision time r is infinite a.s. 
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Proof. We define U = — J2i<j los(Aj — Aj) on t G [0, r]. Applying Ito formula, using fl3.2p 
and the fact that d (Aj, Xj) = 6ij4g'^{Xi)h'^{Xj)dt we obtain 

^^^y dX - ^A, ^ 1 C^(A A)+^(A„A,) ^ ^ ^^(,) ^ ^^(2) + ^^(3) 

^ A,- -A, 2 (A,-A,02 

where 

g{Xi)h{Xi)dvi - g{Xj)h{Xj)duj 



dM = 2Y^ 

i<j 1 



Ai — A,' 



Xj — Xi 

{g\X,)-g\X,))(h\h)-h\^^)) 



~^ (Aj - Ai)2 

^^(3)^v- 1^ V /^(A^_G(A^\ ^ 
A,- — Aj f-7|, .V Aj — A/c A,- — Afc / 

_ sr~^ G{Xj, Xk){Xk — Xj) — G{Xi, Xk){Xk — Xi) + G{Xi, Xj){Xj — Xi) 

(A,-A.)(A.-A„)(A.-A,) • 

We will show that the finite variation part of U is bounded on any interval [0, t\. Lipschitz 
continuity of b, g"^ and h"^ implies that A^^^ < Kp{p — l)t/2 and A^^^ < K'^p{p — l)t, where 
i^T is a constant appearing in the Lipschitz condition. Observe also that if for every x, y, z 
we set 

z) = [(/(x) - g\z)){h\y) - h\z)) + {g\y) - g\z)){h\x) - h\z))]{y - x), 
then H{x, y, z) = {G{x, y) — G{x, z) — G{y, z) + G{z, z)){y — x) and 

H{x, y, z) + H{y, z, x) - H{x, z, y) = 2{z - y)G{y, z) - 2{z - x)G{x, z) 
+ 2{y - x)G{x, y) + G{x, x){z - y) - G{y, y){z - x) + G{z, z){y - x). 

Using the last equality and the fact that \H{x, y, z)\ < 2K^\{y — x){z — y){z — x)\ we can 
write 2rfy4(3) ^ ^^(4) ^ ^^(5)^ ^^leie < A^^ < K^p{p - l){p - 2)t/Q and 

^4(5) _ \^ G{Xj, Xj){Xk — Xi) — G{Xi, Xi){Xk — Xj) — G{Xk, Xk){Xj — Xi) 
' ".1^, (A,-A.)(A.-A.)(A,-A,) 

— I G{Xj, Xj) — G{Xi, Xi) G{Xk, Xk) — G{Xj, Xj)\ 1 

\ Xj — Xi Afc — Xj / Afc — Aj 



i<j<k 



If G{x,x) = 2g'^{x)h'^{x) is convex then obviously the expression under the last sum and 
is non-positive. When G{x,x) is C^'^, (i.e. \G'{x,x) — G'{y,y)\ < G\x — y\) then it 
is bounded by G and \Ai \ < Gt. 

Since finite-variation part of U is finite whenever t is bounded, applying McKean 
argument, we obtain that U can not explode in finite time with positive probability and 
consequently r = oo a.s. □ 

Remark 1. Note that if p = 2 then the assumptions on g'^h'^ can be dropped since in 
that case dA^^^ = 0. Observe also that the Proposition [1] holds also in the complex case. 
The eigenvalues of the process Xt on Tip verify the system (13. 9p and the proof of the 
Proposition [1] remains valid. 
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3.4. Matrix Yamada-Watanabe theorem. 

Theorem 5. Consider the matrix SDE 1\) on Sp 

dXt = g{Xt)dBMXt) + h{X,)dBjg{Xt) + b{Xt)dt 

where g, h,b : H ^ H and Xq G Sp. Suppose that 

\g{x)h{x) - g{y)h{y))\'' < p{\x - y\), x,yeR, (3.14) 

where p : (0, oo) — )■ (0, oo) is a Borel function such that J^.^, p~^{x)dx = oo, that the 
functions G{x,y) = g'^{x)h'^{y) + g'^{y)h^{x) is locally Lipschitz and strictly positive on 
{x ^ y} and that b is locally Lipschitz. Then, for t < t, the pathwise uniqueness holds 
for the SDE / TO) . 

Remark 2. The hypothesis in Theorem [S] on the strict positivity of G{x,y) off the 
diagonal {x = y} is equivalent to the condition that g and h have no more than one zero 
and their zeros are not common. 

Remark 3. In the matrix SDE (13.11) the functions g and h appear only in the martingale 
part, whereas in the equations (13.21) and (13. 3p they intervene also in the finite variation 
part. That is why a Lipschitz condition on the symmetrized function g'^ © h'^ cannot be 
avoided in a matrix Yamada-Watanabe theorem on Sp. 

Proof. We diagonalize Xq = HoXoHq. We first show that the equations (13.21) and (13. 3p 
have unique strong solutions when Aq = Aq and Hq = Hq. The functions 



bi{Xi, . . . , Ap) = b{Xi) + 



Cij{Xi, . . . , Ap, /ill, hi2, . . . , hpp) — dkjhik ^ ^ \ ~i 

Aj — Ak 

dij{Xi, . . . , Ap, /ill, hi2, . . . , hpp) = —-hij ^ ' 2 

^ fc^, {^k - Xj) 

are locally Lipschitz continuous on D = {0 < Ai < A2 < . . . < Ap} x [—1, 1]'', r = p^ 
Thus, they can be extended from the compact sets 



I^m = {0 < Ai < A2 < . . . < Ap < m, Ai+i -Xi> 1/m} x [-1, 1]' 

Lm 



to bounded Lipschitz continuous functions on R^+''. We will denote by 6™, and d^- 



such extensions for m = 1, 2, . . .. 

We consider the following system of SDE (recall that Pkj = Pjk) 

dXT = g{XT)du, + 6r(A")rft, ^ = 1, . . . 

dhij = ^ c^(A'", H) dPkjit) + rf^(A'", H) dt, l<t,j<p. 

Since \g{x)h{x) — g{y)h{y))\'^ < p{\x — y\) and J^_^_ p{x)~^dx = 00, by Theorem [2] with 
q = ^p{p — 1), we obtain that there exists a unique strong solution of the above-given 
system of SDEs. Using the fact that Dm C Dm+i, linim-*oo -Dm = D and the standard 
procedure we get that there exists a unique strong solution (A^, Hf) of the systems (13. 2p 
and (13. 3p up to the first exit time from the set D. This time is almost surely equal to r, 
the first collision time of the eigenvalues. 

Suppose that X and X are two solutions of (13. ip with Xq = Xq = /iqAo/iq . Consider 
the corresponding eigenvalues and eigenvectors processes (A, H) and (A, H) with Aq = 
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Aq = Aq and Hq = Hq = Hq. We have just proved that A(t) = A(t) and H(t) = H{t) for 
every t > a.s., so X{t) = X{t) a.s. □ 

The Theorem [S] and the Proposition [T] imply the following global strong existence result 
for matrix SDEs on the space Sp. 

Corollary 1. Consider functions g,h : H ^ H. Suppose that b,g'^,h'^ are Lipschitz 
continuous, g'^h'^ is convex or in class C^'^ and that G{x, y) is strictly positive on {x ^ y}. 
Then the matrix SDE 1\) on Sp admits a unique strong solution on [0, oo). 



rfXf^t) = 2^Xr{t)dB,{t) + 2{iy + l)dt + AXl''\t)Y, 



Proof. Recall that if a non- negative function F is Lipschitz continuous then y/ F is 1/2- 
Holder continuous. Observe that if g"^ and /i^ are Lipschitz continuous then g'^h^ is locally 
Lipschitz and gh is 1/2-Holder. Thus fl3.14p is verified and the Theorem [5] applies on [0, r). 
By the Proposition [H r = oo almost surely. □ 

4. Applications 

4.1. Noncolliding particle systems of squared Bessel processes. In a recent paper 
by Katori,Tanemura [13], particle systems of squared Bessel processes BESQ^''^ z/ > —1, 
interacting with each other by long ranged repulsive forces are studied. If there are 
particles, their positions xj^'^^ are given by the following system of SDEs, see |14] p. 593: 

dt 

xr\t) - xf\t) 

I x'^^Ht) + x'^'^Ht) 

= 2^J xt\t)dB,it) + 2{u + N)dt + ' \' [ dt, t = l,...,N 

with a collection of independent standard Brownian motions {Bi{t),i = 1, . . . , A^} and, 
if —1 < u < 0, with a reflection wall at the origin. Theorem H] implies that the pro- 
cesses Xl"^ (t) are the eigenvalues of a complex Wishart (or Laguerre) process, with shape 
parameter 6 = u + N, see the end of the Section 3.2. 

Theorem 6. The system of SDEs for a particle system of N squared Bessel processes 
BESQ^''\ with < X[''\'^) < ^(0) < . . . < xi^\o) admits a unique strong solution on 
[0, oo) for u > —1. 

Proof. Like for a Squared Bessel process on R"*", one must start with the following system 
of SDEs 



dYt\t) = 2^\Y}^\t)m{t) + 2{u + N)dt + 2 / i = i,...,n, 

^ Yr\t)-Yl''\t) 

which is well defined on up to the first collision time r. We suppose that < 
y}'"\0) < Yi''\0) < ... < F^''^(0). It follows from the Proposition [U that the collision 
time for the processes (y/'^^(t)), z = l,...,A^isr = oo a.s. 

First suppose that u > —1. The Theorem [2] applied to the last system, with a standard 
use of localization techniques as in the proof of Theorem |5l gives the existence of a 
pathwise unique strong solution {Y-^''\t)). It remains to show that Y}'^\t) > for all 
t > 0. 

Denote 

\Yl''\t)\ + \Yj^\t)\ 



b,{t) = u + N + J2 



Y}''\t)-Yl'\t) 
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We define two stopping times 

= inf{t > 0| Y}''\t) < 0}; 
K = mi{t > bi{t) = 0}. 

Suppose that P('i9 < oo) > 0. Then there exists T > such that P('i9 < T) > 0. As 
Y}''\'d) = and bii'd) = u + N - (N - I) = u + I > 0, we see that if < oo then k > 13. 

It follows from |20|Lemma 3.3, p. 389 that the local time L^(y}'^^) = 0. Using Tanaka's 
formula [20]VI(1.2) we obtain for t > 



J'dAT ^ ^ \ I- s 

= -2E / 1 M 6i(s)rf.<0. 



In the last inequality we used the fact that hi{s) > when < s < k. Thus 

Y}''\{^ + t) AkAT) > 
for t > which contradicts the definition of ^. We deduce that ^ = oo almost surely. 

In the case u = —1, let To = inf{t > 0| Y^''\t) = 0}. Observe that if To < oo, 
then 6i(To) = 0. Define Y}''\t) = Y}''\t) when t < Tq and Y}''\t) = when t > Tq. 
Then {Yl''\Yi''\ . . . , F^'^^) is a solution of the same SDE system as {yI''\ Yi''\ . . . , Y^''^). 
Consequently, by Theorem [2], we have Y^'^^ = Y-^"'' > 0. □ 



4.2. Wishart stochastic differential equations. Wishart processes on Sp are matrix 
analogues of Squared Bessel processes on R"*". Wishart processes with shape parameter 
n (which corresponds to the dimension of a BESQ on R^) are simply constructed as 
Xt = Nf Nt where A'^t is an n x p Brownian matrix. Let a > and B = {Bt)t>o be 
a Brownian p x p matrix. The Wishart stochastic differential equation for a Wishart 
process with a shape parameter a is 

dXt = ^/TtdBt + dBjJXl + aldt 

* (4.1) 

Xq = Xq. 

It was introduced by Bru [3] by first writing the SDE for Xt = Nf Nt and next replacing 
the parameter n by a. It was shown in |3] that if Xq G and a G (p — l,p + 1) then 
there exists a unique weak solution of (14. ip and the condition a > p + 1 implies that 
(14. ip has a unique strong solution. We reinforce considerably these results. 

Our methods apply to the following matrix stochastic differential equation 

dYt = \/\Yt\dBt + dBf^/\Y^\ + aldt (4.2) 

where a & H, Yq = yo & Sp and \Yt\ is defined by taking absolute values of eigenvalues of 
Yt, see the beginning of Section 3. We have g(x) = ^/\x\, h(x) = 1 and G(x, y) = \x\ + \y\ 
for x, 1/ G R. These functions satisfy the hypotheses of the Theorem [5] and the Proposition 

m 
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By Theorem [21 the eigenvalues of the generahzed Wishart process Yf verify the following 
system of SDEs 

dX, = 2^\d,^, + ( « + E ^^^X^ ) 

V kj^i * J 

First, using the Theorem [3] and the Proposition [T] we obtain 

Corollary 2. For a G R and < Ai(0) < A2(0) < . . . < Ap(0) the eigenvalues Aj(t) 
never collide, i.e. the first collision time r = oo almost surely. 

Next, the Theorem [S] implies 

Corollary 3. The generalized Wishart SDE ^4-^ with Yq = uq E Sp has a unique strong 
solution on [0, oo) for any a G R. 

In order to consider the equation fl4.ip . we must prove the non-negativity of the smallest 
eigenvalue of the process F^, when starting from a non-negative value. 

Proposition 2. If a > p — 1 and Ai(0) > then the process \i{t) remains non-negative. 

Proof. We argue as in the proof of the Theorem El □ 

Consequently, using the unicity of solutions in Theorem [5l we obtain 

Corollary 4. The Wishart SDE with xq G S:^ has a unique strong solution on 

[0, oo) for a>p — 1. 

Remark 4. Bru [3] showed that for a > p — 1 the Wishart processes have the absolutely 
continuous Wishart laws which are very important in multivariate statistics, see e.g. the 
monograph of Muirhead [18]. The singular Wishart processes corresponding to a = 
1, . . .p — 1 are obtained as Xt = Nj Nt where A't is an a x p Brownian matrix. Then 
Xq = Nq No has eigenvalue of multiplicity p — a so xq ^ Sp. 

Remark 5. Let Q G S^ . As observed in Bru [3] p. 746- 747, Xt is a solution to the matrix 
SDE fl3.ll) if and only if Yf = QXtQ is a solution of 

dYt = VYtdBtQ + QdBj^/Vt + aQ'^dt. (4.3) 

The process Yt is a Wishart process with shape parameter a and scale parameter . 
Observe that if Q 7^ a/, the matrix SDE (14.31) has not the form (13. ip for any 6, h : R"^ — 
R+. However, using the relation Yt = QXtQ and Theorem |5l it is easy to see that the 
SDE (14. 3p has a unique strong solution when a > p — 1. 

Another perturbation of the Wishart SDE (14. ip is the equation for the Wishart process 
with constant drift c > 0, which may be also viewed as a squared matrix Ornstein- 
Uhlenbeck process 

dXt = ^/XtdBt + dB'[^/x't + aldt + cXtdt, Xq G Sp. (4.4) 

This equation has the form (13. ip with g{x) = a/x, h{x) = 1 and b{x) = ex. By Theorem 
[5l and Proposition [H it has a unique strong solution with t G [0, 00) for any a > p — 1 
and c > 0. More general squared matrix Ornstein-Uhlenbeck processes were first studied 
by Bru [3] and recently by Mayerhofer et al. ^7\. Our strong existence and uniqueness 
result for (14. 4p is not covered by these papers. 
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4.3. Matrix Jacobi processes. Let Op and Ip be zero and identity p x p matrices. 
Define Sp[0, 1] = {X e Sp\ Op < X < Ip}. Denote by ^[0, 1] = {X e Sp\ Op < X < Ip} 
and by 5p[0,/] the set of matrices in 5p[0,/] witli distinct eigenvalues. A matrix Jacobi 
process of dimensions [q, r), witli q Ar > p — 1, and witli values in Sp[0, 1], was defined 
and studied by Doumerc (TUI as a solution of the following matrix SDE, with respect to 
a p X p Brownian matrix Bt 



dXt = yOQdBt,/T~^t+ ^/T^^tdB^VX't + iqI^-{q + r)Xt)dt ^ ^ 

(4.5) 

Xo = xoeSp[0,I]. 

In [To] Th.9.3.1, p. 135 it was shown that i{qAr>p + l and Xq G Sp[0, 1] then (14. 5 P has 
a unique strong solution in Sp[0, 1]. In the case gorr G {p — l,p + 1) and Xq G Sp[0, 1] 
the existence of a unique solution in law was proved in [TU]. Our methods allow one to 
strengthen the results of Doumerc. 

Corollary 5. Let q Ar > p — 1 and xq G iSp[0,/]. Then the matrix SDE has a 

unique strong solution for t G [0, oo). 

Proof. We apply Theorem [5] and Proposition [T] with g[x) = \/\x\, h{x) = a/|1 — x| and 
b{x) = q — {q + r)x. Next we prove similarly as in the proof of the Theorem [6] that 
< Ai(t) < . . . < Ap(t) < 1. □ 

4.4. /3-Wishart and /3-Jacobi processes. Let (3 > 0. One calls a /3-Wishart process a 
solution of the system of SDEs 



dXi 



kdu, + /3 L + ^ ) ^^-^^ 

The /3-Wishart processes were studied by Demni [7j. In the theory of random matrices 
and its physical applications, the /3-Wishart processes are related to Chiral Gaussian En- 
sembles, which were introduced as effective (approximation) theoretical models describing 
energy spectra of quantum particle systems in high energy physics. Usually a symmetry 
of Hamiltonian is imposed and it fixes the value of /3 to be 1, 2 or 4, respectively in real 
symmetric, hermitian and symplectic cases. On the other hand, from the point of view 
of statistical physics, /3 is regarded as the inverse temperature, /3 = l/^ksT), and should 
be treated as a continuous positive parameter. In this sense, the /3-Wishart systems are 
statistical mechanics models of "log-gases" (The strength of the force between particles 
is proportional to the inverse of distances. Then the potential, which is obtained by 
integrating the force, is logarithmic function of the distance. So the system is called a 
"log-gas"). For more information on log-gases, see the recent monograph of Forrester 

In [7] the existence and uniqueness of strong solutions of the SDE system (14.61) was 
established for /3>0,a>p— 1 + and t G [0, r A Rq), where Rq = inf {t| Ai = 0}. Our 
Theorem [2] and Proposition [H together with comparison techniques like in the proof of 
the Theorem El imply the following 

Corollary 6. The SDE system (g7^ with < Ai(0) < A2(0) < . . . < Ap(0) has a unique 
strong solution for t G [0, oo), for any a > p — 1 and /3 > 0. 

The /3- Jacobi processes (Aj),z = 1,. . . ,p are [0, l]^-valued processes generalizing pro- 
cesses of eigenvalues of matrix Jacobi processes defined by (14. 5p : 

dX, = 2y/K{l-K)du, + p(q-{q + r)A, + ^ - >^k) + \k{l - K) \ ^^^^^ 
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Indeed, for /3 = 1 the formula (14. 7p was shown in jTU] and it follows directly from the 
Theorem[31 /3-Jacobi processes were recently studied by Demni in He showed that the 
system (14. 7p has a unique strong solution for all time t when /3 > and qAr > p — 1 + 
As an application of Theorem |2l Proposition [T] and the comparison techniques like in the 
proof of the Theorem [6l we strengthen essentially this result. 

Corollary 7. The SDE system with < Ai(0) < A2(0) < . . . < Aj,(0) < I has a 

unique strong solution for t G [0, oo), for any /3 > and q A r > p — 1. 

Remark 6. It would be interesting to extend our generalization of the Yamada theorem 
to the SDE's considered by Cepa-Lepingle On the other hand the Jacobi eigenvalues 
processes being an important example of the radial Cherednik processes, we conjecture 
that the strong existence and unicity would hold for radial Cherednik processes. For 
radial Dunkl processes this is proved by Demni [6], using j^. 

Acknowledgements. We thank N. Demni, C. Donati, M. Katori and M. Yor for 
discussions and bibliographical indications. 
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